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Abstract 

A  thermodynamic  analysis  is  carried  out  for  electrokinetic  energy  conversion.  We  demonstrate  that  the  efficiencies  depend  solely  on  the 
figure  of  merit  Z  and  are  independent  of  the  working  mode  (generator  or  pump)  at  the  conditions  of  maximum  power  output  and  maximum 
efficiency.  Between  these  two  extreme  points,  the  ratio  of  output  powers  and  the  ratio  of  operating  conditions  (generation  voltage  or  pumping 
pressure)  are  also  functions  of  only  Z  and  independent  of  the  working  mode.  The  figure  of  merit  Z  associated  with  phenomenological 
coefficients  is  less  than  one  due  to  the  intrinsic  entropy  production  in  electrokinetic  flows.  We  establish  the  phenomenological  coefficients 
through  electro-hydrodynamics,  and  find  that  Z  is  dependent  on  the  non-dimensional  electrokinetic  radius,  normalized  zeta  potential  and  non- 
dimensional  parameter  ft  characterizing  the  liquid  property.  The  narrow  range  of  realistic  ft  limits  the  magnitude  of  Z  to  between  0.168  and 
0.465  for  aqueous  solutions.  The  corresponding  maximum  efficiency  thus  varies  between  4.59  and  15.5 1%.  Within  this  range  the  electrokinetic 
devices’  performance  at  maximum  efficiency  is,  however,  close  to  that  at  maximum  power  output.  This  feature  is  identified  in  a  characteristic 
performance  curve  of  efficiency  against  output  power. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

An  electrokinetic  device  can  be  operated  in  either  a  gen¬ 
erator  or  a  pump  mode.  As  a  generator,  a  pressure  difference 
is  imposed  on  a  liquid  in  a  fine  capillary,  and  an  induced 
streaming  current  and  an  induced  electrical  potential  differ¬ 
ence  (called  the  streaming  potential)  are  generated  between 
the  ends  of  the  capillary,  indicating  a  conversion  of  mechan¬ 
ical  energy  into  electrical  energy  (named  as  the  generation 
power).  As  a  pump,  an  electrical  field  is  applied  across  the 
liquid  in  a  capillary  and  if  the  resultant  electro-osmotic  flow 
is  obstructed,  a  pressure  difference  appears  indicating  a  con¬ 
version  of  electrical  energy  into  mechanical  energy  (named 
as  the  pumping  power).  The  streaming  potential  and  electro- 
osmotic  effects  were  first  observed  by  Reuss  [1]  in  1809  and 
Quincke  [2]  in  1859,  respectively.  Both  effects  are  associated 
with  electrical  double  layers  forming  adjacent  to  solid-liquid 
interfaces,  which  were  first  proposed  by  Helmholtz  [3]  in 


*  Corresponding  author.  Tel.:  +1  416  978  1282;  fax:  +1  416  978  7753. 
E-mail  address:  dli@mie.utoronto.ca  (D.  Li). 

0378-7753/$  -  see  front  matter  ©  2005  Elsevier  B.V.  All  rights  reserved, 
doi:  10. 1016/j.jpowsour.2005. 05.057 


1879.  After  that,  several  scientists  such  as  Stern  and  Debye 
have  made  great  contributions  to  the  double-layer  theory  and 
the  electrokinetic  phenomena  (including  electrophoresis). 
A  nice  overview  of  this  development  was  presented  by 
Hunter  [4]. 

The  pioneering  study  of  electrokinetic  energy  conversion 
might  owe  to  Osterle  [5]  who  predicted  in  1964  a  conversion 
efficiency  of  0.392%  with  water  as  the  working  liquid.  Since 
then,  dozens  of  papers  devoted  to  electrokinetic  energy 
converters  have  been  published.  Electrokinetic  pumps  have 
received  most  of  the  interest  because  they  can  generate 
high  pressures  and/or  high  flow  rates  for  chip-based  fluid 
manipulation  like  microelectronics  cooling  [6]  and  chro¬ 
matographic  separation  [7].  These  pumps  are  either  based 
on  glass  columns  packed  with  silica  microspheres  [8-14] 
or  composed  of  parallel  “slot”  structures  micro-fabricated 
in  glass  [15]  and  silicon  [16]  substrates.  Santiago  and  co¬ 
workers  [17,18]  surveyed  the  progress  in  the  development 
of  microscale  electrokinetic  pumps.  The  experimentally 
measured  maximum  pumping  efficiency  is  mostly  less  than 
1%,  and  seems  much  lower  than  the  theoretically  predicted 
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Nomenclature 

a 

capillary  radius  (m) 

e 

charge  of  proton  (1.602  x  10-19  C) 

E 

applied  electrical  field  (V  m-1) 

8i 

defined  non-dimensional  coefficients,  i=  1, 

Q 

2, 

G 

C) 

hydrodynamic  conductance  (m4  s-1  Pa-1) 

i 

electrical  current  density  (Am-2) 

I 

electric  current  (A) 

Boltzmann’s  constant  (1.3807  x  10-23  JK~ 

6 

K 

non-dimensional  electrokinetic  radius 

l 

capillary  length  (m) 

no 

bulk  ionic  density  (m-3) 

M 

phenomenological  coefficient  characterizing 

the  electro-osmotic  flow  (m4  V-1  s-1) 

P 

hydrodynamic  pressure  (Pa) 

P 

applied  pressure  gradient  (Pam-1) 

Q 

liquid  flow  rate  (m3  s-1) 

r 

radial  coordinate  (m) 

R 

universal  gas  constant  (8.314  J  mol-1  K-1) 

S 

electrical  conductance  (Sm) 

T 

liquid  temperature  (K) 

u 

liquid  electrokinetic  flow  velocity  (ms-1) 

W 

output  power  (W) 

Z 

figure  of  merit 

Zy 

valence  of  ions 

Greek  letters 

p 

non-dimensional  property  of  liquid 

1 1 

conversion  efficiency 

£ 

dielectric  constant  of  liquid  (C  V-1  m-1) 

0 

electrical  potential  (V) 

A 

molar  conductivity  of  liquid  (m2  S  mol-1) 

ti 

liquid  viscosity  (Kgm-1  s-1) 

Pe 

net  charge  density  (C  m-3) 

(70 

electrical  conductivity  of  liquid  (S  m-1) 

\jf 

electrical  double-layer  potential  (V) 

$ 

zeta  potential  of  the  capillary  surface  (V) 

£ 

non-dimensional  zeta  potential  of  the  capillary 

surface 

Subscripts 

max  W 

maximum  output  power 

max  rj 

maximum  conversion  efficiency 

efficiency.  Without  restricting  to  small  surface  potentials 
or  thin  double-layers,  numerically  calculated  pumping 
efficiencies  could  be  10%  [19]  and  15%  [20]  in  a  cylindrical 
capillary  at  the  conditions  of  maximum  pumping  power  and 
maximum  pumping  efficiency,  respectively. 

Only  a  few  of  past  papers  have  investigated  electroki- 
netic  generators.  Forty  years  ago  Burgreen  and  Nakache 


[21]  predicted  an  efficiency  of  approximately  17%  for 
a  generator  using  water  in  equilibrium  with  the  CO2  in 
the  air.  However,  the  maximum  efficiency  of  a  recently 
developed  electrokinetic  microchannel  battery  using  tap 
water  is  barely  on  the  order  of  0.01%  [22].  More  recently, 
Daiguji  et  al.  [23]  investigated  the  electro-chemo-mechanical 
energy  conversion  in  nanofluidic  channels  and  predicted  an 
efficiency  of  more  than  1%.  Another  discovery  is  a  voltage 
induced  by  a  liquid  flow  in  carbon  nanotube  bundles  [24]. 
While  the  underlying  mechanism  is  still  unclear  [25,26],  this 
phenomenon  explores  the  potential  of  energy  conversion  in 
nanotubes. 

In  analyzing  the  efficiency  of  electrokinetic  energy 
conversion,  the  great  majority  of  previous  studies  employed 
directly  the  theory  of  electrokinetic  microchannel  flow. 
As  numerical  solutions  were  necessarily  involved  in  these 
analyses,  however,  the  true  factors  affecting  the  devices  per¬ 
formance  were  missed  or  left  behind  the  complex  formulae. 
In  essence,  electrokinetic  devices  utilize  the  liquid  flow  and 
electrical  current  to  realize  the  thermodynamic  conversion 
between  mechanical  energy  and  electrical  energy.  Therefore, 
we  turn  to  a  thermodynamic  analysis  of  electrokinetic  energy 
conversion  in  this  paper,  and  provide  a  new  understanding  of 
devices  performance  at  the  conditions  of  maximum  power 
output  and  maximum  efficiency.  We  noticed  that  Osterle  and 
co-workers  had  employed  the  phenomenological  equations 
to  investigate  the  electrokinetic  energy  conversion  [27,28]. 
However,  they  focused  mainly  on  the  maximum  conversion 
efficiency,  and  ignored  the  performance  at  the  maximum 
power  output.  Therefore,  the  comparison  of  device  per¬ 
formance  and  operating  conditions  at  the  two  extreme 
points  is  still  lacking  from  the  thermodynamic  viewpoint. 
Moreover,  the  parametric  study  of  conversion  efficiency 
and  output  power  is  incomplete  from  the  electrokinetic 
viewpoint. 

This  paper  is  organized  as  follows.  In  Section  2  a 
thermodynamic  analysis  is  presented  for  irreversible 
electrokinetic  energy  conversion  process.  The  operating 
conditions  and  devices  performance  at  maximum  power 
output  and  maximum  efficiency  are  derived  in  terms  of 
phenomenological  coefficients,  and  an  important  parameter 
Z,  named  as  figure  of  merit,  is  defined.  In  Section  3  the  phe¬ 
nomenological  coefficients  and  also  Z  are  specified  through 
electro-hydrodynamics.  We  examine  in  Section  4  the  factors 
affecting  Z  and  thus  the  performance  of  a  practical  electroki¬ 
netic  generator  as  an  example.  A  summary  of  our  results  and 
a  few  technical  comments  are  provided  in  Section  5. 

2.  Thermodynamic  analysis 

Within  the  linear  response  regime  of  electrokinetic  flow, 
the  Onsager  reciprocal  theorem  relates  the  liquid  flow  rate 
Q  and  the  electrical  current  I  though  a  capillary  (of  any 
cross-sectional  shape),  to  the  pressure  difference  A p  and  the 
electrical  potential  difference  A 0  between  its  entrance  and 
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exit  [29,30]: 

Q  =  G(—Ap)  +  M(—Acp)  (1) 

/  =  M(-Ap)  +  S(-A0)  (2) 

where  G  represents  the  hydrodynamic  conductance,  M  char¬ 
acterizes  the  electro-osmotic  flow  in  Eq.  (1)  and  the  streaming 
current  in  Eq.  (2),  and  S  indicates  the  electrical  conduc¬ 
tance.  It  is  noted  that  in  the  generation  mode,  the  pressure 
decreases  in  the  direction  of  fluid  flow  (to  provide  mechani¬ 
cal  power)  and  the  voltage  increases  in  the  direction  of  current 
(acting  as  a  battery).  In  the  pumping  mode,  however,  the  pres¬ 
sure  increases  in  the  direction  of  fluid  flow  and  the  voltage 
decreases  in  the  direction  of  current  (to  provide  electrical 
power).  As  such,  the  output  power  W  and  the  conversion  effi¬ 
ciency  r)  are,  respectively,  given  by: 

W  =  IAcp  for  generators,  and  W  =  QAp  for  pumps 

(3) 

IAcp  QAp 

T)  =  - for  generators,  and  r\  =  - for  pumps. 

Q(-Ap)  I(-A<f> ) 

(4) 

Substituting  Q  and  I  in  Eqs.  (3)  and  (4)  with  the  expressions 
in  Eqs.  (1)  and  (2)  yields  W  and  ij  that  are  both  functions 
of  solely  the  generation  voltage  A <p  for  generators  (where  a 
known  pressure  difference  Ap  is  applied))  or  the  pumping 
pressure  Ap  for  pumps  (where  a  known  potential  difference 
Acp  is  applied).  Therefore,  we  can  maximize  W  and  ij  sep¬ 
arately  by  differentiating  them  with  respect  to  either  A0  or 
Ap.  This  mathematical  manipulation  is  straightforward  and 
hence  omitted  here.  The  resultant  performance  formulae  are 
summarized  below. 


2.7.  Maximum  generation  power  or  pumping  power 


For  electrokinetic  generators,  the  generation  power  is 
maximized  when  the  generation  voltage  achieves: 

M 

A0maxW  —  ^v;(— Ap) 

where  the  subscript  max  W  signifies  the  condition  of  max¬ 
imum  generation  power  (or  the  maximum  pumping  power 
for  liquid  pumps  presented  later).  Note  that  A cpmSLXw  in  Eq. 
(5)  is  exactly  half  of  the  maximum  voltage  produced  by  the 
generator,  when  the  electrical  current  I  in  Eq.  (2)  vanishes 
and  hence  W  and  ij  are  both  zero.  The  maximum  generation 
power  and  the  corresponding  conversion  efficiency  are  given 
by: 


fkmax  W  —  ^ZG(A p) 


7 max  W  — 


z 


2(2  -  Z) 


M2 
GS ' 


(6) 

(7) 

(8) 


One  can  see  that  ijmaxw  is  solely  dependent  on  the  non- 
dimensional  parameter  Z  that  is  less  than  1  according  to 
Prigogine’s  theorem  [31].  Within  this  limit,  ?7max  ^  is  a  mono- 
tonically  increasing  function  of  Z  and  we  shall  call  it  the 
“figure  of  merit”  following  Osterle  and  co-workers  [27,28]. 
In  fact,  this  figure  of  merit  Z  can  be  viewed  as  the  product  of 
the  streaming  current  effect  in  Eq.  (9)  and  the  electro- osmotic 
effect  in  Eq.  (10)  (refer  to  Eqs.  (1)  and  (2))  [29,31]: 


(9) 

(10) 


These  two  effects  occurring  simultaneously  in  electrokinetic 
devices  underlie  the  energy  conversion  principle  of  electroki¬ 
netic  generators  and  pumps.  The  significance  of  the  figure  of 
merit  Z  will  be  explored  in  the  next  section. 

In  the  pumping  mode,  the  maximum  pumping  power  and 
the  corresponding  pumping  pressure  are  expressed  as: 


WmaxW  =  lzS(A0)2 


(11) 


A/?max  W 


M 

2G 


(—Acp). 


(12) 


As  this  condition,  however,  the  pumping  efficiency  is  identi¬ 
cal  to  the  generation  efficiency  in  Eq.  (7). 


2.2.  Maximum  generation  efficiency  or  pumping 
efficiency 


For  electrokinetic  generators,  the  conversion  efficiency  is 
maximized  when  the  generation  voltage  satisfies: 

A</>max t]  =  ^(1  -  Vl  -  Z)(— A0)  (13) 

where  the  subscript  max  r]  indicates  the  condition  of  maxi¬ 
mum  generation  or  pumping  efficiency.  At  this  condition,  the 
generation  power  and  efficiency  are  given  by: 

vTzzz(i  -  vTzzz)2  , 

wmaxv  =  G- - - - -(Ap)2  (14) 


_  d- 

7 max  rj  — 

Eq.  (15)  also  applies  to  the  maximum  pumping  efficiency 
when  the  pumping  pressure  reaches: 

Apmuxtj  =  — (1  -  Vl  -  Z)(— A<t>).  (16) 

The  corresponding  pumping  power  is  described  by: 

vTz^za  -  vnvz)2  7 

wmax„  =  S- - - LiM)  2.  (17) 

Osterle  [5]  first  proved  that  the  efficiency  of  electrokinetic 
energy  converters  is  independent  of  the  working  mode,  which 


VT^zy 


(15) 
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is  attributed  to  the  Onsager  reciprocal  relation  (see  Eqs.  (1) 
and  (2)).  Evidently  the  maximum  generation  or  pumping  effi¬ 
ciency  in  Eq.  (15)  is  also  a  monotonically  increasing  function 
of  the  figure  of  merit  Z.  Moreover,  we  notice  that  both  the 
ratio  of  output  powers,  Eq.  (18),  and  the  ratio  of  operating 
conditions,  Eq.  (19),  between  the  two  extreme  points  dis¬ 
cussed  above  are  also  functions  of  solely  Z  and  independent 
of  the  working  mode: 

wmaxw  (l  +  ^/^^z)2 

-  =  - ,  -  (to) 

ffmax  rj  4>/ 1  —  Z 


A/?max  W  _  A0max  \y  _  1  +  \/I  Z  (19) 

APmax  rj  A 0m ax  r)  2 

As  seen  from  these  two  equations,  the  figure  of  merit  Z  gauges 
the  performance  of  electrokinetic  energy  converters.  The  ulti¬ 
mate  goal  of  electrokinetic  devices  researchers  is  to  maximize 
Z  in  Eq.  (8),  which  can  be  realized  by  improving  the  stream¬ 
ing  current  effect  in  Eq.  (9)  and/or  the  electro-osmotic  effect 
in  Eq.  (10).  This  issue  will  be  addressed  later  in  this  paper. 

Fig.  1  compares  the  efficiencies  under  the  conditions  of 
maximum  output  power  and  maximum  efficiency.  As  stated 
previously,  both  efficiencies  increase  with  the  rise  of  Z.  How¬ 
ever,  only  at  Z  >  0.6  can  these  two  efficiencies  be  distin¬ 
guished  (5%  difference).  In  the  high  Z  limit  (i.e.,  Z^  1), 
OmaxW  achieves  0.5  while  Omaxri  approaches  1.  We  will 
demonstrate  later  in  this  paper  that  the  realistic  Z  is  outside 
of  this  region.  Therefore,  the  above  two  extreme  conditions 
are  essentially  close  to  each  other  in  typical  electrokinetic 
energy  converters.  The  two  ratios  defined  in  Eqs.  (18)  and 
(19)  are  demonstrated  in  Fig.  2  against  different  values  of  Z. 
The  discrepancy  in  the  output  powers  at  the  two  extremes  is 
again  less  than  5%  unless  Z>0.6.  However,  the  operating 
conditions,  say  generation  voltage,  are  significantly  different 
provided  that  Z  >  0.2.  The  magnitude  of  generation  voltage 
(or  pumping  pressure,  not  shown  in  the  axis  label  of  right  ordi¬ 
nate)  at  the  condition  of  maximum  output  power  is  always 


Figure  of  merit,  Z 


Fig.  1.  Efficiencies  vs.  the  figure  of  merit  Z  under  conditions  of  maximum 
output  power  and  maximum  conversion  efficiency  for  electrokinetic  energy 
converters. 
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Fig.  2.  Ratios  of  output  power  and  generation  voltage  for  generators  (or 
pumping  pressure  for  pumps)  vs.  the  figure  of  merit  Z  between  the  conditions 
of  maximum  output  power  and  maximum  conversion  efficiency. 


lower  than  that  at  maximum  efficiency.  In  the  high  Z  limit 
(i.e.,  Z  — >  1 ),  the  former  is  only  half  of  the  latter.  Note  that 
77max  w  is  also  half  of  r/max  r]  in  this  limit. 


3.  Mathematical  formulation 


We  have  demonstrated  in  the  preceding  section  that  the 
electrokinetic  devices  performance  is  associated  with  the 
coefficients  G,  M  and  S  in  the  phenomenological  Eqs.  (1) 
and  (2).  In  this  section,  we  want  to  specify  these  coefficients 
through  the  electrokinetic  flow  theory  for  the  capillary  model, 
i.e.,  a  cylindrical  tube  of  length  /  and  radius  a.  However,  the 
same  method  presented  here  also  applies  to  channels  of  other 
cross-sectional  shapes  such  as  slit  or  rectangular  channels.  It 
is  well  known  that  the  liquid  electrokinetic  flow  velocity  u 
consists  of  electro-osmotic  flow  (the  first  term  on  the  right 
hand  side  of  Eq.  (20))  and  pressure-driven  flow  (the  second 
hand  on  the  right  hand  side  of  Eq.  (20))  [32]: 


sf  f  0  \  A0  a 2  /  r2  \  A p 

/a  \  f  )  l  4/x  \  a2 )  l 


(20) 


where  e  is  dielectric  constant  of  the  liquid,  /x  the  liquid  vis¬ 
cosity,  r  the  radial  coordinate,  0  the  electrical  double-layer 
potential  induced  by  the  surface  charge  at  the  channel  wall 
that  is  often  characterized  by  the  zeta  potential  f .  Neglecting 
the  diffusive  ion  fluxes  in  the  direction  of  the  applied  electri¬ 
cal  field  and  assuming  a  Boltzmann  distribution  of  ions,  the 
local  current  density  i  in  a  symmetric  electrolyte  is  given  by 
[4,19,20]: 


i  =  ao  cosh(zve0/k#r) 


+  PqU 


(21) 


where  ctq  is  the  electrical  conductivity  of  the  bulk  liquid, 
Zy  the  valence  of  ions,  e  the  charge  of  proton,  k#  the  Boltz¬ 
mann’s  constant,  T  the  liquid  temperature,  and  pe  the  net 
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charge  density  that  is  determined  from  the  Poisson  equation: 
Pe  =  -eVV-  (22) 


Integrating  Eqs.  (20)  and  (21)  over  the  channel  cross-section 
and  comparing  with  the  phenomenological  Eqs.  (1)  and  (2) 
yield  [33]: 


4 

7ta 
8  ill 


(23) 


M  =  — 


Tt£a2f 


111 


g  l 


5  = 


7t£2^2 

III 


(g2  +  Pg3) 


(24) 

(25) 


gi 


82 


=  1-2 


d  WQ 
d  (r/a) 


83  = 


(26) 

(27) 

(28) 


It  is  noted  that  gi  defined  in  Eqs.  (26)-(28)  are  all  positive 
and  independent  of  the  sign  of  zeta  potential.  In  these 
equations,  the  distribution  of  double-layer  potential  \j/  is 
solved  from  the  Poisson-Boltzmann  equation  [4]: 


V 


2 


2zYeno 


sinh 


£ 


(29) 


where  no  is  the  bulk  ionic  density.  To  this  end,  the  figure  of 
merit  Z  in  Eq.  (8)  can  be  specified  as: 


ggf 

82  +  A?3  ' 


(30) 


There  appear  three  important  non-dimensional  parameters  in 
the  definition  of  Zin  Eq.  (30).  The  first  one  ft  is  the  ratio  of  the 
conductive  current,  A  in  Eq.  (31),  to  the  convective  current 
due  to  electro-osmosis,  £//i  in  Eq.  (31),  and  essentially  char¬ 
acterizes  a  non-dimensional  property  of  the  working  liquid1 : 


P 


Ail 

Jrt 


(31) 


where  A  is  the  molar  conductivity,  and  R  the  universal  gas 
constant.  As  gi  in  the  definition  of  the  figure  of  merit  Z, 
Eq.  (30),  are  positive,  a  smaller  p  is  desirable  in  order  to 
attain  a  higher  Z.  Since  A  varies  only  between  about  4  and 
20mSm2mol_1  for  aqueous  solutions  [34],  however,  the 
realistic  values  of  ft  likely  span  the  small  range  2</3<  10. 
We  will  demonstrate  in  the  next  section  that  it  is  mainly  this 
small  range  of  /3  that  limits  the  magnitude  of  Z.  The  other 


1  In  Griffiths  and  Nilson’s  paper  [19],  the  reciprocal  of  ft  is  named  as 
Levine  number. 


two  parameters  that  hide  in  the  definitions  of  gi  (see  Eqs. 
(26)-(28))  are  the  non-dimensional  electrokinetic  radius: 


2z2e2fto 

£k^T 


(32) 


and  the  normalized  zeta  potential: 


(33) 


Up  to  this  point,  the  phenomenological  coefficients  and  the 
figure  of  merit  Z  have  been  determined.  The  only  problem 
remaining  to  be  solved  is  the  distribution  of  double-layer 
potential  in  Eq.  (29),  which  will  be  addressed  in  the 
proceeding  section. 


4.  Results  and  discussion 

To  examine  the  effects  of  the  three  non-dimensional 
parameters  p ,  K  and  f  on  electrokinetic  devices  perfor¬ 
mance,  we  first  need  to  determine  the  double-layer  potential 
x/j.  Here,  Eq.  (29)  was  analytically  solved  with  the  method 
proposed  by  Levine  et  al.  [35]  that  is  applicable  to  cylindri¬ 
cal  channels  with  arbitrary  zeta  potential  and  arbitrary  non- 
dimensional  electrokinetic  radius.  The  integrations  involved 
in  Eqs.  (26)-(28)  were  numerically  evaluated  with  Simpson’s 
rule. 

Fig.  3  shows  a  characteristic  performance  curve,  i.e., 
efficiency  r)  against  output  power  W,  of  an  electroki¬ 
netic  generator  working  with  pure  water  (of  10~5M 
ionic  concentration).  The  generation  power  has  been  nor¬ 
malized  by  its  maximum  value  WmSLXw-  The  capillary 
radius  is  a  =  0.2  pan.  Other  necessary  parameters  are  [34]: 


Fig.  3.  Characteristic  performance  curve  of  electrokinetic  generators.  The 
generation  power  is  normalized  by  its  maximum  value.  The  curve  vanishes 
in  the  limits  of  low  (A)  and  high  (D)  generation  voltage  due  to  the  vanishing 
voltage  and  the  vanishing  electrical  current,  respectively.  The  two  adjacent 
points  are  maximum  generation  power  (B)  and  maximum  generation  effi¬ 
ciency  (C).  Pure  water  (of  10-5  M  concentration)  is  the  working  liquid,  and 
the  capillary  radius  is  0.2  pun.  All  other  parameters  are  referred  to  the  text. 
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viscosity  /x  =  0.9  x  10-3  Kgm-1  s_1,  relative  dielectric  con¬ 
stant  s  =  80  x  8.854  x  10“ 12  C  V-1  m_1,  molar  conductivity 
A  =  10  m  S  m2  mol-1 ,  zeta  potential  f  =  — 100  mV,  and  room 
temperature  T-  298  K.  The  resultant  figure  of  merit  Z  is  0.264 
where  the  involved  non-dimensional  parameter  ft  is  5.128. 
The  four  points  noted  on  the  curve  are:  (1)  A  represents  the 
vanishing  efficiency  in  the  low-generation  voltage  limit;  (2) 
B  represents  the  maximum  generation  power;  (3)  C  repre¬ 
sents  the  maximum  generation  efficiency;  (4)  D  represents 
the  vanishing  efficiency  in  the  high-generation  voltage  limit. 
The  points  A  and  D  collapse  to  the  origin  because  either  the 
generation  voltage  (at  point  A)  or  the  electrical  current  (at 
point  D)  is  zero  so  that  the  generation  power  vanishes.  The 
two  extreme  points  B  and  C  are  close  to  each  other  because  the 
figure  of  merit  Z  under  investigation  is  small  (refer  to  Fig.  1). 

We  are  interested  in  the  likely  magnitude  of  Z  when  /3 
is  within  the  practical  range  of  2  <  /3  <  10  for  typically  used 
aqueous  solutions.  Fig.  4  shows  the  contours  of  Z  as  a  func¬ 
tion  of  the  non-dimensional  electrokinetic  radius  K  and  the 
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Fig.  4.  The  contours  of  figure  of  merit  Z  as  a  function  of  the  non-dimensional 
electrokinetic  radius  K  and  the  normalized  zeta  potential  £*  at  ft  =  2  (top) 
and  f>-  10  (bottom).  The  highest  magnitudes  of  Z  are,  respectively,  0.465 
(top)  and  0.168  (bottom). 
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Fig.  5.  The  contours  of  maximum  generation  efficiency  ijmaxr]  as  a  func¬ 
tion  of  the  non-dimensional  electrokinetic  radius  K  and  the  normalized  zeta 
potential  at  ft  =  4.  The  highest  magnitude  of  rjm axr]  is  9.66%. 

normalized  zeta  potential  £  at  /3  =  2  (top)  and  /3=  10  (bot¬ 
tom),  respectively.  In  this  range  of  /3,  the  highest  magnitude 
of  Z  varies  only  between  0.465  (top)  and  0.168  (bottom), 
within  which  the  performances  at  the  conditions  of  maxi¬ 
mum  power  output  and  maximum  efficiency  are  essentially 
indistinguishable  (refer  to  Fig.  1).  Interestingly  and  impor- 
tantly,  Z  achieves  its  summit  in  a  range  of  K  and  f  ,  not  a 
particular  set  of  these  two  parameters.  For  example  at  /3  =  10 
(bottom  subfigure  in  Fig.  4),  Z  remains  nearly  constant  at 
f  >4  despite  that  K  drifts  toward  lower  values  as  f  grows 
larger.  In  other  words,  it  is  unnecessary  to  choose  extremely 
high  zeta  potentials  and/or  extremely  small  pores.  This  pre¬ 
diction  is  consistent  with  those  from  Griffiths  and  Nilson 
[19]  and  Min  et  al.  [20],  who  analyzed  the  performances 
at  maximum  pumping  power  and  maximum  pumping  effi¬ 
ciency,  respectively.  We  also  notice  that  as  decreases,  the 
lower  limit  of  £  that  optimizes  Z  becomes  smaller  while  the 
corresponding  K  remains  almost  unchanged. 

As  the  generation  or  pumping  efficiencies  at  the  two 
extreme  conditions  increase  with  the  rise  of  the  figure  of 
merit  Z  (see  Fig.  1),  they  should  follow  the  same  pattern  as 
Z  increases  with  K  and  f*.  Fig.  5  shows  the  contour  of  max¬ 
imum  conversion  efficiency  at  /3  =  4,  which  is  quite  similar 
to  the  efficiency  contour  at  the  condition  of  maximum  out¬ 
put  power  (not  shown).  By  choosing  optimum  values  of  K 
and  £*,  we  can  get  an  efficiency  of  9.66%  at  most.  If  /3  =  2  is 
assumed,  instead,  the  maximum  efficiency  could  go  as  high 
as  15.51%.  This  efficiency  in  close  agreement  with  Min  et 
al.  [20]  analysis  is  much  higher  than  the  best  experimen¬ 
tal  efficiency  reported  to  date  [18].  It  turns  out  that  /3  is  the 
determinant  factor  in  the  electrokinetic  conversion  efficiency. 
Therefore,  more  attention  should  be  paid  to  the  search  of 
electrolyte  solution  with  a  high  /3  in  the  future.  For  example, 
Kirby  and  co-workers  [36,37]  have  recently  demonstrated  the 
performance  improvement  of  electrokinetic  micropumps  by 
using  zwitterionic  solute  additives. 
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Other  than  the  figure  of  merit  Z,  the  output  powers  under 
the  conditions  of  maximum  power  output  and  maximum  con¬ 
version  efficiency  are  also  related  to  the  liquid  conductance 
(specifically,  hydrodynamic  conductance  for  generators  and 
electrical  conductance  for  pumps,  see  Eqs.  (6),  (1 1),  (14)  and 
(17)).  However,  these  two  powers  are  approximately  the  same 
in  realistic  electrokinetic  devices  due  to  the  small  magnitude 
of  Z.  Hence,  we  investigate  only  the  maximum  generation 
power  in  Eq.  (6)  and  the  maximum  pumping  power  in  Eq. 
(11)  that  are  specified  as: 


ffma x  W 
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4 

2  na 

ZP -  for  generators,  and 

32  ji 

2ns2  E2  2  2 
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Fig.  6.  The  contours  of  maximum  generation  power  (top)  and  maximum 
pumping  power  (bottom)  as  a  function  of  the  non-dimensional  electrokinetic 
radius  K  and  the  normalized  zeta  potential  at  f  =  4.  The  applied  pressure 
gradient  is  P=  1  atm  cm-1  for  the  generator  and  the  applied  electrical  field 
is  E  =  100  V  cm- 1  for  the  pump.  The  unit  for  the  contour  bar  is  nW  m-1  for 
both  cases. 


where  P=—Ap/l  and  E-—  A  0/ /  denote  the  applied  pressure 
gradient  and  electrical  field,  respectively.  One  can  see  that  the 
maximum  generation  power  is  a  strong  function  of  capillary 
radius  and  thus  a  larger  pore  is  desirable  as  well  as  a  higher 
figure  of  merit  Z.  The  maximum  pumping  power  is,  however, 
not  directly  dependent  on  Z  and  thus  not  limited  by  the  non- 
dimensional  liquid  property  ft.  Yao  and  Santiago  [17]  have 
demonstrated  numerically  that  g i  (see  Eq.  (26))  increases 
with  the  rise  of  non-dimensional  electrokinetic  radius  K. 
Also,  g i  is  only  a  weak  function  of  zeta  potential  f  (or  f '). 
Therefore,  the  maximum  pumping  power  is  enhanced  when 
K  and/or  £  become  larger.  Fig.  6  shows  the  contours  of  max¬ 
imum  generation  power  (top)  and  maximum  pumping  power 
(bottom)  at  /3  =  4.  The  applied  pressure  gradient  for  the  gener¬ 
ator  is  P  =  1  atm  cm- 1 ,  and  the  applied  electrical  field  for  the 
pump  is  E-  100  V  cm-1 .  As  expected,  both  powers  increase 
with  the  rise  of  non-dimensional  electrokinetic  radius  K  (i.e., 
capillary  radius).  As  the  maximum  generation  power  (top)  is 
proportional  to  Z,  it  follows  the  similar  trend  to  Fig.  4  against 
the  normalized  zeta  potential  f  .  For  liquid  pumps,  however, 
the  maximum  pumping  power  is  elevated  at  higher  f  . 


5.  Conclusions 

Based  on  the  general  Onsager  relation  in  electrokinetic 
flow,  we  have  conducted  a  thermodynamic  analysis  for 
electrokinetic  generators  and  pumps.  Their  conversion 
efficiencies  at  the  condition  of  maximum  power  output  and 
the  condition  of  maximum  efficiency  are  exactly  the  same. 
These  two  efficiencies  depend  solely  on  the  figure  of  merit  Z. 
This  Z  also  governs  the  ratio  of  output  powers  and  the  ratio 
of  operating  conditions  (i.e.,  generation  voltage  or  pumping 
pressure)  between  the  two  extreme  points.  Moreover,  these 
two  ratios  are  independent  of  the  working  mode.  The  figure 
of  merit  Z  expressed  by  the  phenomenological  coefficients 
is  less  than  1  due  to  the  intrinsic  entropy  production  in 
irreversible  electrokinetic  flows.  In  the  high  limit  of  Z,  the 
maximum  efficiency  approaches  1  with  zero  power  output 
while  the  efficiency  at  maximum  power  output  achieves 
0.5.  In  the  same  limit,  the  generation  voltage  or  pumping 
pressure  at  maximum  power  output  is  reduced  to  half  of  that 
at  maximum  efficiency. 

We  have  also  determined  the  phenomenological 
coefficients  and  thus  the  figure  of  merit  Z  through  electro¬ 
hydrodynamic  analysis.  We  find  that  Z  is  dependent  on 
non-dimensional  electrokinetic  radius  K ,  normalized 
zeta  potential  f  and  non-dimensional  liquid  property 

This  /3,  the  ratio  of  conductive  current  to  convective 
current  due  to  electro-osmosis,  spans  the  small  range 
2</3<10  for  aqueous  solutions.  As  a  result,  the  optimal 
Z  calculated  numerically  varies  only  between  0.465  and 
0.168.  The  corresponding  maximum  efficiency  varies 
between  4.59  and  15.51%.  The  latter  value  is  much  higher 
than  the  best  experimentally  reported  efficiency  of  5.6% 
[36],  indicating  a  potentially  large  developing  space  for 
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state-of-the-art  electrokinetic  pumps.  We  also  find  that  the 
practically  small  Z  has  dimmed  the  performance  difference 
at  maximum  power  output  and  maximum  efficiency.  This 
feature  has  been  demonstrated  in  the  performance  curve 
of  efficiency  against  power  output  for  an  electrokinetic 
generator  working  with  pure  water.  In  addition,  our  analysis 
suggests  that  future  work  on  electrokinetic  energy  convert¬ 
ers  should  be  focused  on  the  search  of  electrolytes  with 
a  high  /3. 
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